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NOTATION 


Free-surface  amplitude  at  a  point  (x,y,0)  per  unit  amplitude  of 
oscillation  in  the  kth  mode 


Reflected  wave  amplitude  per  unit  incident  wave  amplitude 


Free-surface  amplitude  per  unit  incident  wave  amplitude  at  a 
point  on  the  incident  side  of  the  wave  due  to  Incident  and 
reflected  waves 


Free-surface  amplitude  per  unit  incident  wave  amplitude  at  a 
point  on  the  lee  side  of  the  wave  due  to  transmitted  wave 


A(I) 


Free-surface  amplitude  at  a  point  on  the  Incident  side  per  unit 
incident  wave  amplitude  due  to  Incident,  reflected  and  motion- 
generated  waves 


Free-surface  amplitude  per  unit  Incident  wave  amplitude  at  a 
point  (x,y,0)  due  to  combined  sway,  heave,  and  roll  motion 


'ki 


Added  mass  in  the  kth  mode  due  to  motion  in  the  fcth  mode 


b  Half-beam  of  a  cross  section  at  the  calm  waterline 

Damping  (due  to  wave  making)  coefficient  in  the  kth  mode  due  to 
motion  in  the  ith  mode 


Restoring  coefficient  (spring  constant)  in  the  kth  mode  due  to 
motion  in  the  fcth  mode 


F.(e)  Wave-excited  force  or  moment  in  the  kth  mode 

k 

g  Gravitational  acceleration 

J  -  Imaginary  unit 

K  -  A  /  2tt  Wave  number 

N  Two-dimensional  unit  normal  vector 

n  Unit  normal  vector  on  the  body  surface  pointing  into  the  body 

0,x,y,z  Right-handed  Cartesian  coordinate  system  (see  Figure  1) 

U  Forward  velocity  of  ship 


v 


C(x,y,t) 

«A 

X 

V 

CR(x,y,z,t) 

£y(x,y,z,t) 

P 

<f(x,y,z,t) 

4>D(x,y,z) 

^(x.y.z) 

4>k(x,y,z) 

<f»k(y.z,x) 

<Ps(x,y,z) 


Free-surface  elevation 
Incident  wave  amplitude 

Have  length 

Have  heading  angle;  y  •  0  for  following  wave 

Complex  amplitude  of  displacement  of  body  due  to  oscillation 
in  the  kth  mode 

Vertical  displacement  of  a  point  (xty,z)  on  a  body  relative  to 
the  point  on  the  free  surface  on  the  same  vertical  line 
(see  Equation  (6)) 

Vertical  displacement  of  a  point  (x»y,z)  on  a  body  due  to  motion 


Density  of  water 

Velocity  potential  function  which  represents  the  total  fluid 
disturbance  due  to  wave  and  body  motion 

Complex  velocity  potential  for  diffracted  wave 
Complex  velocity  potential  for  incident  wave 

Complex  velocity  potential  for  forced  oscillation  in  the  kth  mode 
Two-dimensional  complex  velocity  potential 


Velocity  potential  for  disturbance  generated  by  the  steady 
forward  motion  of  a  ship 
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Abstract 

An  analytical  method  Is  developed  for  pre¬ 
dicting  the  vertical  motion  of  a  point  on  a 
ship  relative  to  the  motion  of  the  free  sur¬ 
face.  The  prediction  method  presented  here 
takes  into  account  the  effect  of  the  deforma¬ 
tion  of  the  incident  waves  on  the  relative 
motion.  The  causes  of  the  deformation  consid¬ 
ered  are  the  waves  generated  by  diffraction  and 
the  waves  generated  by  the  motion  of  the  ship. 
The  method  is  based  on  the  two-dimensional 
approximation  of  the  flow  around  the  cross 
sections  of  ships.  The  results  reveal  that  the 
deformation  of  incident  waves  Is  so  significant 
that  it  should  be  accounted  for  in  the  pre¬ 
diction  of  the  relative  laotlon  of  ships. 

I.  Introduction 

The  vertical  motion  of  a  point  on  a  ship 
hull  with  respect  to  the  undulating  free  sur¬ 
face  is  Important  Information  In  the -seakeeping 
Investigation  of  ships.  This  motion  Is  often 
called  "Relative  Motion."  The  relative  nation 
has  a  direct  effect  on  the  Inception  of  deck 
wetness,  slamming  of  the  ship  bottom,  and 
rudder  and  propeller  emergence . 

In  general,  the  relative  motion  is  computed 
under  the  assumption  that  the  Incident  wave 
system  is  undisturbed.  However,  the  incident 
waves  can  be  significantly  disturbed  In  the 
vicinity  of  a  ship  due  to  the  diffraction  by 
the  ship  surface  and  the  waves  generated  by  the 
motion  of  the  ship.  Hence,  one  can  easily 
surmise  that  the  cause  of  the  poor  correlation 
between  the  predicted  and  the  measured  '»^,3 
relative  motion  Is  the  assunption  of  the  un¬ 
disturbed  incident  waves  near  a  ship. 

In  this  paper  a  method  to  account  for  the 
free-surfece  disturbance  In  the  computation  of 
relative  motion  Is  described.  The  method,  as 
an  initial  attempt.  Is  limited  to  a  two-dimen¬ 
sional  approximation  within  the  context  of 
strip  theory.  This  approach  Is  taken  because 
firstly,  the  results  can  be  readily  Incorpo¬ 
rated  into  the  existing  computational  scheme  of 
ship  motion  based  on  strip  theory*  and  secondly, 
an  evaluation  of  the  two-dimensional  approxi¬ 
mation  ought  to  be  made  before  undertaking  a 


more  complex  three-dimensional  approach.  The 
strip  approach  of  obtaining  the  free-surface 
disturbances  near  a  ship  hull  was  encouraged 
by  the  success  achieved  by  strip  theory  In 
the  computation  of  the  absolute  motion  of  ships 
in  waves. 

The  two-dimensional  potentials  are  obtained 
by  using  the  method  of  distribution  of  pulsat¬ 
ing  sources  on  the  boundary  of  the  cross 
section  of  the  body.  The  source  distribution 
is  extended  on  the  waterline  inside  the  body  to 
remove  the  Irregular  behaviors  of  the  poten¬ 
tials  at  certain  discrete  freguencies.  Various 
cross  checkings  of  the  numerical  convergence 
are  made  to  ensure  the  validity  of  the  computed 
results. 

The  computed  results  of  a  pontoon  having  a 
uniform  cross  section  identical  to  the  midship 
of  a  mariner  hull  form  are  presented.  The 
contribution  from  the  various  sources  generat¬ 
ing  local  waves  near  the  pontoon  are  shown  In 
the  figures.  The  results  indicate  clearly  that 
the  free  surface  motion  at  the  sides  of  a  ship 
should  be  considered  If  a  reliable  prediction 
of  the  relative  vertical  motion  is  desired. 

The  present  study  will  be  incorporated  in  a 
strip  fashion  into  an  existing  ship-motion 
computer  program,  and  the  validity  of  the  pre¬ 
sently  developed  method  will  be  investigated. 

II.  Analysis 

The  coordinate  system  to  be  used  in  the 
analysis  Is  a  right-handed  Cartesian  coordinate 
system  which  moves  on  the  calm-water  surface  in 
the  mean  course  of  the  ship  with  the  ship  speed. 
The  origin  Is  located  directly  above  the  center 
of  gravity  of  the  «hlp  at  its  mean  position; 
the  x-axis  is  directed  toward  the  bow;  and  the 
z-axls  it  directed  vertically  upward;  see 
Figure  I. 
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We  assuae  that  the  eater  I*  Incompressible 
and  Its  motion  Irrotatlonal  such  that  a  velo¬ 
city  potential  can  be  defined  In  the  fluid 
region.  We  alto  attune  that  the  water  It  in¬ 
finitely  deep  and  that  no  current  exists. 

If  we  denote  the  velocity  potential  which 
represents  the  disturbance  of  the  fluid  by 
0(x,y,z,t),  then  It  can  be  decoaiposed  In  the 
following  form  for  a  ship  moving  with  a 
constant  velocity  II  In  a  regular  plane  wave 
system. 


♦(x.y.z.t)  -  -Ux  ♦  (x.y.z) 

Jw  t 

♦  He(B,(x,y,z)e  *  1  (1) 

Here,  Is  the  wave-encounter  frequency,  Ha 
means  the  real  part  of  what  follows,  f,  Is  the 
steady  potential,  and  f,  is  the  oscillatory 
complex  amplitude  of  tne  oscillatory  potential 
which  can  bn  further  decomposed  Into 

6 

♦.  •  *i  ♦  ♦£  i2) 

where  ft  Is  the  incident-wave  potential;  fg 
the  dlffracted-wave  potential;  Sk  the  forced- 
oscillation  natentlal  In  the  kth  mode  of 
motion;  and  (fc  the  «x*plex  amplitude  of  the 
displacement  of  the  body  due  to  oscillation  In 
the  kth  mode.  The  incident-weve  potential  dj 
can  be  given  explicitly  by 

♦jlx.y.i,  -  .-Jk(xcosu  ♦  yslnw)  ♦  Kz  (3) 
where 


"a 

*  J  f 


e  •  •;,*  c,e 


J<*,t 


M 


where 


C#(x,y)  -  »°)  ^ 

c,(*.y)  ■  •j-  ^ 

The  vertical  displacement  of  a  point  (x,y,z) 
on  a  ship,  denoted  by  Ey(x,t),  Is  given  by 

Ev(*,t)  -  *  ♦  E,(t)  -  xE,(t)  ♦  yEjt)  (5) 


where  k  •  (x,y,z),  and  Ej.  E»  and  E.  are, 
respecTlvely,  the  heave,  roll  and  pitch  dis¬ 
placement  from  the  naan  position  of  the  body. 
The  relative  motion  of  a  point  with  respect  to 
the  free  surface  motion  at  the  same  horizontal 
coordinates  (x,y)  Is  defined  by 

E*(x.t)  -  Ev(x,t)  -  C(x.y.t)  (6) 

The  vertical  position  of  the  point  on  the 
hull  from  the  calm  water  surface  at  zero  speed 
is  not,  of  course,  necessarily  the  same  as  that 
at  a  non-zero  forward  speed  due  to  the  sinkage 
and  trim  of  the  ship.  In  a  strict  sense,  an 
Inclusion  of  the  sinkage  and  trim  effect  in  the 
determination  of  the  relative  motion  means  that 
the  second-order  effects  contributed  by  the 
terms  such  as  0(psp().  0(d|) ,  and  0(pj)  should 
also  be  Included  In  the  analysis.  However,  an 
investigation  of  the  second-order  effects  will 
not  be  pursued  In  the  present  study.  The 
analysis  will  be  kept  within  the  first  order 
of  the  Incident  wave  amplitude  and  of  the 
slenderness  parameter  of  the  body. 


u  •  wave  frequency  In  radians  par  sec 
g  •  gravitational  acceleration 
CA  *  wave  amplitude 

K  •  w*/g  •  2w/X  -  wave  maker  for  deep 
water 


Determination  of  Steady  Free  Surface 

The  free-surface  deformation  caused  by  a 
steady  translation  of  a  ship  in  calm  water  at 
constant  speed  U  is  obtained  from  the  line¬ 
arized  Bernoulli's  equation  as 


X  -  wave  length 

v  -  wave  heading  angle;  |i  •  0  lq  fallowing 
waves 

J  -  /=T 


Cs(x,y)  -  Jj-  P^x.y.O)  (7) 

The  boundary- value  problem  for  ♦  is  as 
follows: 


The  free-surface  elevation  c(x,y,t)  can  be 
obtained  from  Bernoulli's  equation  In  terms  of 
the  velocity  potential  by 


t(*.y.t) 


-  7  {ft  *  ♦  •<•*> 

g-  *  ♦sx(",r'0)*jM*t] 


When  the  spatial  variables  x,  y ,  «  and  n 
1 normal)  are  used  as  a  subscript.  It  means  a 
partial  dtrlvatlso  with  the  respective  veri- 
able.  Also,  when  a  complex  function  Is 
multiplied  by  o'***,  only  the  Heal  part  of 
the  product  should  bo  real  lead. 


♦sx*U'v'0)  *  $  m  0  (8) 


where  S,  denotes  the  hull  surface  below  the 
calm-weter  surface,  and  n  •  (n. ,nt,n.)  is  the 
unit  normal  vector  on  S(~polntlng  Into  the  body. 


and 

♦#f(«,y.-«)  ■  o 

(10) 

1  im  A  • 

i  Of r~'/2)  for  r  •  0 

ahara 

loll)  for  x  «  0 
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To  obtain  the  solution  of  the  foregoing 
boundary-value  problem,  the  thln-shlp  assump¬ 
tion  will  be  used.  If  we  represent  the  hull 
geometry  by 

y-*f(x,z).  (12) 

then  the  unit  normal  vector  on  the  ship  hull 
surface  pointing  into  the  body  can  be  deter¬ 
mined  by 


<f  .tMJ 

n  •  — —  1  ■  —  -  on  y  ■  *f  (13) 

-  sr+r+rr 

X  z 

Substitution  of  Equation  03)  into  Equation 
(9)  yields 

z‘uf*  (,<,) 

We  assume  that  the  ship  in  consideration  is 
thin  such  that  f  •  0(e)  where  e  is  a  small 
parameter  representing  the  beam  to  length  ratio 
of  the  ship.  Then,  from  Equation  (14),  dis¬ 
carding  the  terms  of  0(c*),  we  find  that 


ct(*.v)  •  ~  ♦„(|t»v,o) 


(x-t)f.(C.C) 

- - rtr  «dC 

[(*-U,ey*>C*]VI 


Uf  (,6) 

s<»> 

„  If  we  substitute  the  expression  for 
Gox(*~£>y>0  from  Equation  (17)  into  Equation 
( 18) ,  we  get 


Cs(x.y) 


fe(e.C)d?dc 


r 


dQsecG 


.  /••slnik^lcose]  cos(ky  sin6)kek;dk 
•'o  k  *  $  secl0 


♦1Y(*.±0.z)  -  TOf^  (15) 

The  solution  of  ♦,  Is  well  known  from 
chin-ship  theory  (see,  e.g.,  Mehausen  and 
Lai tone$)  and  is  given  by 

U  ft  M5.0 

♦.(n.y.z)  -  2=1/  »  -  : . .  d£dC 

2v^  Ax-ow^(*-«)* 


-  Tiff  G,(*-e.V.*K>f5(e,C)dW<;  06) 

S<») 

where  denotes  the  longitudinal  center- 

plane  of  the  ship,  and 


G0(x-C,y,z*c)  -  -  ===== 


0  0 

.  coslk(x-C)cos81cos(kv  sin9) 
kcos’0  -  g/U1 


2W  ,  A(z+Oseca0  _ 

*  iHI  d6sec*9  eu  sln(jj(x-5)sec8) 


3/ 


cos  1^1  ysinSsac'e)  (17) 


where^ means  the  principal -value  Integral. 
Thus,  from  Equations  (7)  and  (18),  we  get 


2*  jJ|C*ecle 

dOsec30  eu 


&ff  k  6lAf 

s(.)  0 

■c°s[^(x-E)sece|  co*(£  y  s  in*0secJ  ()j  (19) 


Within  the  first-order  approximation,  the 
wave  profile  along  the  side  of  the  hull  can  be 
obtained  by 


dOsecB 


tj(x.o)  ■  «el-  U//  fcdr.dc/* 1,1 
1  «(•)  0 


*/  dk  &// 


/2s 


•I  desec’e  ekosec 


*ez  J 


(20) 


where  ky  -  4,  Z’  •  t  ♦  i  (x-E)cos6,  and 
{  ■  ,  or 


4S(*,0)  “  *•  I*  JT ff  f^dC 

S<*> 

•  |  f  d6sec fl/*  ek2  dk 

■'o  •'o 

/2*  «■  . 
*«c9d9/  -%,-,ecSe  dk} 
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*  faff  V*dt/?'d8sec,e«k,”cJ0Z'J  (2I) 

«!•>  o  j 


which  amounts  to  the  wave  profile  along  the 
longitudinal  centerplane. 

Determination  of  Unsteady  Free  Surface 


The  free-surface  deformation  caused  by  the 
incident  waves,  diffracted  waves,  and  mot  Ion- 
generated  waves  Is  obtained  from  Equations  (kb) 
and  (2).  The  unknown  functions  are  ^  and 
f .  In  the  following  sections  we  describe  pro¬ 
cedures  for  obtaining  these  unknown  quantities. 

Diffraction  Potential.  From  the  kinematic 
boundary  condition  on  the  hull  surface,  we 
obtain  from  Equation  (3) 

♦onl.  •-♦ml, 

-K(jn,cosu  +  jn2sinp  -  njJgJ  (22) 


If  we  assume  that  the  ship  is  slender  such 
that  n,<<n, ,  n,,  then,  discarding  n,  in 
Equation  (22),  we  obtain 

♦onl  *  K(jn2Sinp  -  n,)^ 


-  j^IjnjSlnu  -  n,) 

.  e-jh(x0cosw  +  y0sinu)  ♦  Kz,  (23) 

where  (x0,y,,z0)  indicates  a  point  on  the  hull 
surface  S&. 

From  Equation  (23)  we  can  Infer  that 

♦D(*.y.z)  =  *(y,z;x)e'jK*cos‘1  (2k) 

where  x  affects  <p  as  a  parameter  rather  than 
an  independent  variable.  By  applying  the 
Laplacian  operator  to  the  right-hand  side  of 
Equation  (2k),  we  obtain 

♦yy  *  *zz  '  K'cos'ud  -  0  (25) 

An  appropriate  linearized  free-surface 
condition  for  the  velocity  potential  6(y,z;x) 
i  s 

(j“e  *  UJx)  *(T'°;*)« -JKxcosu 

♦  «♦,  .'J*XC°*’i  -  0  (26) 

or  using  the  relation  •  u  -  KUcosu,  we  have 


♦zlv.O;x)  -  K6  «  0  (27) 

The  kinematic  body-boundary  condition  for  ♦  is 

♦J  -  juk.(jH2sinp  -  N, ) 

NIC,(x)  A 

.  e-jKyoS'hW  ♦  Kz„  (28) 


where  II  ■  (M2  ,N 
vector  on  C,(x) 
of  a  cross  sect 
t i on  as  |y|  ■* • 


j)  is  the  two-dimensional  unit 
which  is  the  imnersed  contour 
ion  at  x.  The  far  field  condi- 
is  given  as' 


j«JK|y| 
~  1 1  y  I 


In  fact,  we 
problem  for  d  i 
hence ,  I f  we  f i 
function,  G,  wh 
(27)  and  (29),  i 
the  form  of 

♦(y.t;x)  - 


si ny  .  , 

1  for  u  d  x 

(29 

for  u  *  n 

find  that  the  boundary-value 
s  confined  in  the  y-z  plane; 
nd  an  appropriate  Green's 
ilch  satisfies  Equations  (25), 
we  can  obtain  the  solution  in 


/« 


I ) G ( I  ;y,z)dl 


where  the  unknown  function  Q  should  be  found 
from  the  remaining  boundary  condition  (28). 

The  appropriate  Green's  function  was  given  by 
Ursell,!j  and  the  solution  for  ♦  was  obtained 
by  Choo^  and  Troesch.8 

The  approximation  of  by  Equation  (2k), 
where  the  variable  x  is  suppressed  as  a  para¬ 
meter,  has  led  to  the  free  surface  condition 
(27)  and  the  far-field  condition  (29).  As  one 
can  readily  observe,  the  forward-speed  effect 
is  nonexistent  in  the  foregoing  problem.  The 
foregoing  simplification  can  be  criticized  for 
the  lack  of  consistency  in  the  perturbation 
scheme.  However,  the  analysis  will  be  pursued 
on  this  basis  with  the  assumption  that  the 
forward  speed  has  no  significant  effect  on  the 
wave  diffraction.  Furthermore,  under  the 
assumption  that  the  wave  lengths  of  interest 
are  greater  than  the  order  of  ship  beam,  the 
Helndtolz  equation  given  by  (25)  will  be  re¬ 
placed  by  the  Laplace  Equation,  i.e., 

•  +  ♦  -  0  (31) 

yy  zz 

which  is  exact  for  p  ■  tr/2,  i.e.,  beam  waves. 

A  further  digression  from  the  original  problem 
will  be  made  with  a  heuristic  argument,  based 
on  the  success  of  strip  theory,  that  the  near- 
fleld  solution  derived  with  the  radiation 
condition 

llm  (♦  *  JKd)  -  0  (32) 

y*  ’ 

Is  acceptable.* 

t - 

Although  the  diffraction  potential  is  not 
.explicitly  solved  In  the  usual  ship  iwtion 
theory  based  on  the  strip  assumption,  a  sim¬ 
ilar  radiation  condition  is  Invoked  in  apply¬ 
ing  the  Hesklnd  relation  In  the  two  dimen¬ 
sional  sense. 
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Thus,  the  modified  boundary-vein*  problem 
now  is  reduced  to  er.  almost  Identical  problem 
of  forced  oscillation  of  two-dimensional 
cylinders  in  a  free  surface  in  the  combined 
modes  of  heave  and  sway.  It  should  ba  noted 
that  the  diffraction  potential  Is  a  function 
of  the  wave  frequency  u.  This  Is  a  noticeable 
difference  from  the  strip  solution  used  In  the 
ship  motion  computation*  In  which  the  diffrac¬ 
tion  potential  is  treated  as  a  function  of 
wave-encounter  frequency  i%.  In  Appendix  I, 
the  process  of  solving  the  diffraction 
potential  satisfying  the  two-dimensional 
Laplace  Equation  ()l)  with  the  boundary  condi¬ 
tions  (27),  (28)  and  (32)  Is  described.  The 
description  is  given  in  a  general  form  for 
any  prescribed  function  replacing  the  right- 
hand  side  of  Equation  (28). 

notion  and  Radiation  Potential.  The  mot  Ion 
of  the  ship  Is  obtained  by  solving  two  sets  of 
linearized  coupled  equations  of  motion,  which, 
according  to  the  coordinate  system  given  In 
Figure  I.  for  k,t  •  1,  3,  $  and  k.t  ■  2,  *,  8, 
are 

+  ♦«««»-  Fk'’  <”> 

l. 


where  Re.  and  lm,  are  the  real  and  imaginary 
parts,  respectively,  of  a  complex  function, 
the  I mag  I  nary  part  of  which  Is  preceded  by 
j  •  /-I  ;/Ldx  Is  the  integral  owpr  the  length 
In  the  positive  x  direction  and/,  .dl  the 
Integral  over  the  submerged  contOoP'of  the 
cross  section  located  at  x;  b/is  the  two- 
dimensional  approximation  of  g.  ;  and  6,.  is 
the  Kronecker  delta. 

The  restoring  coefficients  C.£  are  given  by 
Cm  ’  P*f.  C,*  -  -pg*.  C,,  -  pg*w,  C,s  -  C,, 

"  P9\,.  Cj,  -  pgW Ixx^f  -08).  »nd  CS5  -  pglf 
* (l yy /If  -  08)  where  V  is  the  displaced  volume. 

Am  the  weterp I ene  area,  M*  the  moment  of  the 
weterplane  area  about  the  y-axis,  lxx  and  lyy, 
respectively,  the  mass  moment  of  inertia  about 
the  x-  and  y-axis,  and  OB  the  vertical  distance 
from  the  center  of  buoyancy  to  the  calm  water¬ 
line.  The  wave  excited  forces  are  given  by 

Fk  ^ 

*'l  "t(x)  1  e 

♦A.  -  Vk,}&]*i <*> 


Here, 

M(mass  of  ship)  for  kwt^3 

I kk (mass  moment  of  inertia)  for 
k-l>  b 

0  for  kft  except  for  H,,  • 

*  -  *,H  where  z,  Is  the 
vertical  coordinate  of  the 
center  of  gravity; 

Akt  is  the  added  mass  coefficient  In  the  kth 
mode  due  to  the  motion  In  the  Eth  mode; 
the  Corresponding  damping  coefficient;  C-,  the 
restoring  coefficient;  and  FJ»)  the  wave* 
excited  force  or  moment.  * 

The  coefficients  In  the  equations  of  motion 
are  obtained  by  slander-body  strip  theory.* 

The  expressions  are  given  as  follows: 

\l  *"*)['.£ /*/  <4kd(y 


for  k  ■  1 ,  2 . 6,  tdiere  h,  »  yNi  -  zN* , 

N,  “  -xN,  and  N(  •  xN,. 

From  the  foregoing  equations  it  is  obvious 
that  If  the  motion  potentials  bk(y,z;x)  are 
known,  we  can  solve  the  equations  of  motion 
and  find  the  six-degrees  of  freedom  motion  £ 
The  first  three  of  Ek  are  the  linear  displace¬ 
ments  from  the  mean  position  of  the  ship  in  the 
x,  y,  and  z  directions,  respectively,  and  the 
remaining  three  are  the  angular  displacements 
about  the  x,  y.  and  z  axes,  respectively.  In 
niaxerlcal  order  of  k  these  are  called  surge, 
swey,  heave,  roll,  pitch  and  yaw  motion. 
Utilizing  the  derivation  shown  in  Appendix  I, 
we  can  show  that 

6  ' 

♦  bD(y,*;x)  b.  (y,z;x)£ 

k-l 

•♦j  ♦  *0* 


(**)  Substitution  of  Equation  (37)  Into  Equation 

*  1  i  (4b)  yields 


skf  *  '-j  [*£/d*/  (*mi  *  2w,As 

1  *  l  C(x) 

♦  **,«*>♦;  -•]  < 


for  k,E  "  1,  2,  ....  6. 


C,(x.y)  ■  -  (bj  (x.Y.O)  ♦  b0) 

?  <♦&  -  ♦’,?.>  - J-?  ♦  «T.) 

.  -  *,>)  (s« 
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III.  Numerical  Checks 

For  the  purpose  of  illustrating  computa¬ 
tions,  a  pontoon  having  uniform  cross  sections 
of  the  midship  of  a  Mariner  Class  ship  is 
chosen.  The  length  of  the  pontoon  Is  arbi¬ 
trarily  taken  as  five  times  the  beam.  The 
cross  section  Is  close  to  a  rectangle  except 
for  the  rounded  bilges.  The  Immersed  divi¬ 
sion  of  the  cross  section  Is  23.06  m  beam  and 
9.07  m  draft.  The  offsets  used  to  describe 
the  section  are  given  in  Table  1  and  the 
section  view  Is  given  In  Figure  2. 

Table  1  -  Offsets  of  Right-Half  of  Midship  Section 
of  Mariner  Claes  Ship 


V 

Z 

(ml 

0. 

-9.000 

4  363 

-9.001 

•706 

-9.00B 

10.300 

-S.40B 

*0.968 

-7.040 

11  400 

-8.629 

11.534 

-4.191 

11.53 

0. 

Figure  2  -  Right-Half  Cross  Section  View  of 
Midship  of  Merinsr  Ctaos  Ship 


As  described  In  Appendix  I,  the  velocity 
potentials  are  obtained  by  using  the  method  of 
source  distribution  on  the  imsersed  contour  of 
the  cross  section  of  the  body.  As  shown  by 
Equation  (1*7),  an  approximation  Is  made  by 
assuming  a  constant  source  strength  on  each 
line  segment  which  makes  up  the  contour.  The 
line  segments  are  made  by  connecting  the  adja¬ 
cent  two  points  on  the  contour  by  a  straight 
line.  The  points  chosen  are  as  shown  In 
Figure  2.  It  can  be  expected  that  numerical 
accuracy  will  increase  as  the  number  of  points 
chosen  on  the  contour  Increases;  however,  a 
compromise  should  be  made  to  minimize  computer 
costs.  A  total  of  eight  boundary  points  are 
taken  on  the  right  half  of  the  Imsersed  contour 
of  a  cross  section  at  shown  by  the  black  dots 
In  Figure  2.  The  boundary  points  thus  chosen 
provided  the  desired  accuracy  end  yielded  a 
satisfactory  numerical  convergence. 


It  is  well  known  that  the  method  of  Green's 
function  to  solve  for  the  velocity  potential 
associated  with  an  oscillating  body  in  a  free 
surface  suffers  from  the  existence  of  indef¬ 
inite  solutions  at  certain  frequencies. 9 
Descriptions  of  the  existence  of  the  Indef¬ 
inite  solutions  and  their  removal  are  given  in 
Appendix  II.  To  remove  the  Indefinite  solu¬ 
tions,  the  source  distribution  Is  extended  on 
the  line  z  •  0  inside  the  body,  and  a  rigid 
wall  condition,  d*  «  0,  Is  Imposed  on  that 
line,  which  will  be  referred  to  as  "top  deck" 
hereafter.  According  to  Frank'O  the  "Irregular 
frequencies"  for  a  rectangle  of  beam  B  and 
draft  T  are  obtained  by 

Kmb  “  “mb^9  ”  T  eoth(m*T/B),  m  *  I,  2,  3,  .  .  . 

where  b  “  y.  For  the  section  considered  here, 
•/T  “  2.56,  and  therefore  we  get  K(b  *  1.86, 
and  K.b  ■  3-19.  In  Figure  3,  the  singular 
behavior  of  the  amplitude  of  the  free  surface 
motion  due  to  beam  regular  waves  at  y  •  1.05b 
(the  wave  Incident  side)  at  the  first  irregular 
frequency  1s  shown.  The  amplitude  of  the  wave 
Is  normalized  by  the  incident-wave  amplitude, 
!•••»  I C p I where  C(  is  obtained  by  Equation 
(38).  The  solid  line  is  obtained  by  taking  two 
line  segments  on  the  top  deck  between  y  »  0  and 
y  *  b,  and  ^  *  0  and  y  »  -b,  and  imposing  the 
condition  $z  -  0.  In  Table  2,  a  comparison  of 
the  free-surface  amplitudes  is  presented  for 
various  top-deck  conditions.  The  values  for 
the  rigid-wall  conditions  imposed  on  up  to  six 
line  segments  on  the  top  deck  are  given  to¬ 
gether  with  that  obtained  by  imposing  the 
♦  “  0  condition  on  the  two  line  segments  on 

the  top  deck.  Also  shown  are  the  values  ob¬ 
tained  by  ■  4  on  the  two  line  segments  on 
the  top  deck  at  Kb  »  0.3,  1.0  and  2.0.  The 
latter  values  are  shown  to  check  Otvnatsu's 
statement'1  that  any  arbitrary  values  for  pj 
on  the  top  deck,  provided  an  appropriate 
syimsetric  or  ant  I -symmetric  condition  Is  main¬ 
tained,  can  be  chosen  to  remove  the  irregular 
frequencies.  It  appears  that  the  nmthod  holds 
In  higher  frequencies  only  in  the  present  case. 
It  looksobvious  from  Equation  (II-6)  that, 
unless  fi  ■  0  or  ♦'  ■  0  Is  imposed  on  the  top 
deck,  there  Is  no  guarantee  that  the  trivial 
solution  for  g’  in  the  Interior  domain  will  be 
obtained.  One  would  assume  that  more  segments 
on  the  top  deck  should  yield  better  results; 
however,  the  values  obtained  near  Kb  ■  1.0  for 
*  and  6  segnmnts  on  the  top  deck  stww  that  the 
assumption  Is  not  necessarily  true.  Why  more 
than  two  segments  on  the  top  deck  show  an 
Irregular  behavior  In  the  vicinity  of  Kb  •  1.0 
Is  not  yet  clear.  Based  on  this  investigation, 
the  two  segments  on  the  top  deck,  together  with 
the  other  segments  on  the  body  contour,  are 
chosen  for  the  ensuing  computations. 
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Table  2  -  >m  Eurtaoa  Amptaudas  My*  1«k 
tar  Variout  Toa-Oach  CondMone 


Kb 

No  Top 

26ag 

6 

♦  *  ® 
28eg 

*1  -  4 

2  Sag 

0.3 

0932 

0.932 

0.932 

0.933 

0.926 

0.642 

os 

1463 

1.428 

1.318 

1.360 

1.610 

0.9 

1528 

1.470 

1.381 

1366 

1.608 

1.0 

1.212 

1.188 

0.680 

0.706 

1.247 

0.848 

1.1 

1.313 

1.271 

2.962 

3034 

1.324 

1.2 

1.946 

1.677 

2.387 

2  427 

1.838 

1.6 

2.143 

1.982 

2.010 

2.009 

2.083 

166 

2168 

1.986 

1.989 

1.986 

1.986 

1.70 

2.179 

1 .977 

1.970 

1.964 

1.988 

1.75 

2.228 

1.970 

1.963 

1.944 

1.982 

1.80 

2  366 

1.962 

1.937 

1.926 

1.976 

186 

3.386 

1.966 

1.922 

1.909 

1.988 

1.90 

1.231 

1.948 

1.906 

1.894 

1.961 

1  96 

1.734 

1.941 

1.897 

1.880 

1.964 

2.0 

1.843 

1  936 

1.887 

1.888 

1.947 

1.937 

2.6 

1.986 

1.877 

1.863 

1.830 

1.883 

Tabta  3  —  Compartaon  of  Wave  Excited  fmaa 
Obtained  by  Eq.OM  and  Eq.  09) 


IF'J'I/A 

If*?  1  /* 

if*:' 

I/A 

Kb 

(361 

1381 

1381 

(381 

(36) 

(30) 

0  1 

01525 

0.1536 

0.8121 

0.8175 

0  0062 

0.0062 

05 

0.9861 

0.5888 

0.5116 

05161 

00602 

0.0500 

0.9 

05710 

0.5761 

0.3429 

0.3449 

00670 

00868 

1.4 

0  4696 

0.4776 

0.2177 

0.2101 

00690 

00692 

20 

03780 

0.3903 

0.1069 

01104 

0  0642 

00854 

3  0  b 


2.6 
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Kb 


Figure  3  Rmovd  af  Sin^pder 
Irregular  Frequency 


Several  cross  checks  of  the  numerical  re¬ 
sults  were  made  to  ensure  the  accuracy  of  the 
computations.  One  was  a  correlation  of  the 
wave-excited  forces  computed  by  two  independent 
methods;  one  by  a  direct  computation  using  the 
expression  given  by  Equation  (36)  and  the  other 
by  an  indirect  method  using  the  damping  coef¬ 
ficients.  The  latter  method  is  given  by 
Newman 12  as 


KiA. 


(pgJ#kk/u) 


t 


(39) 


for  two-dimensional  bodies  for  |i  *  90  degrees. 
The  comparison  of  the  two  independent  results 
given  in  Table  3  shews  a  good  agreement  sdtlch 
gradual ly  deteriorates  up  to  about  3  percent 
difference  at  Kb  -  2.  In  any  case  a  3  percent 
difference  Is  well  within  the  numerical  errors 
resulting  from  the  segmentation  of  the  boundary 
contour  of  a  cross  section.  The  denominator  A 
in  Table  3  equals  to  2pgbtA  for  |FJ*'|and  |F>* 
and  2p9b*cA  For  |F*«>|. 

The  hydrodynamic  coefficients  resulting 
from  the  coupling  between  the  sway  and  roll 
modes  are  Aj»,  A,,,  |tt,  and  From 

Gram's  theorem  we  have 


ff%* In  -  Vkn>*  *  0  (W) 

S 


where  S  Is  the  surface  bounding  the  fluid 
domain  which  consists  of  the  body  surface  S§, 
the  free  surface  Sp,  and  a  vertical  cylinder 
surface  Sr  of  a  large  radius  ft  with  a  bottom 
closure  Sg.  Since  ♦.  «  pk2  ■  ft*,  on  Sr, 

■  0  on  Sr,  and  0(T/*1f)  on  Sr,  we  can 
show  that 

ff  Vitn*  “fj  ♦And* 

»•  S| 

which  means  that 


from  Equations  (34)  and  (3$)-  The  values  of 
these  coefficients  are  shown  In  Table  4, 
Although  the  percent  difference  between  the 
two  values  at  the  higher  frequencies  appears 
large,  the  error  is  within  the  bound  which  is 
associated  with  the  number  of  segments  chosen 
in  the  present  case. 


Table  4  -  Comparison  of  Coupled  Sway-ftoft 
1  AjjSnd  1],l  and  Weft  levy  llg  and 

Coefficient* 


Kb 

*3. 

03 

0.1272 

0  12*6 

0.0934 

0.0316 

1.0 

0  0603 

0.0614 

0.0866 

0.0666 

20 

0  0060 

0  0066 

0.0667 

0.06B 

The  expression  for  a  pulsating  source  of 
unit  strength  below  e  free  surface  Is  given  by 
Equation  (l-5).  The  principal-value  Integral 
in  this  equation  can  be  converted  to  an  expo¬ 
nential  integral  In  the  form 

f  ek<«+<>co»h(y-n)  » 

"'o 

-  «e|-e*loSEI(-l«S)  1  l*e',fllS]; 

the  plus  sign  for  y-ri  >  0  and  the  minus  sign 

for  y-n  <  0 

where 

s  •  (y-n)  ♦  iw 

E.(s)  *r  dt 

■'s 

"  f.i \n«n 

■  T  ♦  In  S  Hrnl  *  l*r#  S!  <  *• 
n*>l 

y  •  O.S772I5W  ...»  Euler's  constant. 


For  the  evaluation  of  E,(S),  one  can  use  the 
above  equation;  however,  when  jsl  is  large, 
the  Infinite  series  makes  the  computation  in¬ 
efficient.  Thus,  one  can  use  various  rational 
and  polynomial  approximations  f or  a  large  argu 
ment.’3  In  this  work  the  Laguerre  quadrature 
method  introduced  by  Todd1*1  is  used. 

To  obtain  the  free  surface  elevation 
f.„(*.v)  the  velocity  potentials  should  be 
evaluated  at  *  •  0,  which  in  turn  means  that 
the  Green's  function  given  by  Equation  (1-5) 
should  be  evaluated  at  z  ■  0.  Then,  for  the 
sources  at  the  top  deck,  the  principal-value 
integral  becomes 

/cosh(y-u)  -  isinkiy-q)^ 
o-k 

•  cosax  CKaX)  ♦  sinaX  |  y  ♦  SI  (aX)| 

-  ijsinaXC(aX)  -  cosaX  |-j-  ♦  SI(aX)|| 


where  X  •  y-n,  and  Si  and  Ci  are  the  sine  and 
cosine  Integrals,  respectively,  which  are  given 
by 


Since  these  two  functions  approach  zero  in  an 
oscillatory  manner  as  x- ”,  the  series  computa¬ 
tion  for  large  x  Is  expected  to  require  a 
large  number  of  terms  before  converging  to  the 
desired  accuracy  (error  less  than  10'’).  The 
rational  approximation  of  Hastings'^  (see, 
e.g.,  p.  233  of  deference  13)  was  employed  to 
obtain  the  sine  and  cosine  integrals,  and  the 
results  were  found  to  be  Identical  to  those  of 
the  laguerre  method  of  Todd. 

The  amplitude  of  the  radiating  waves  at 
y  -  1 00  b  generated  by  the  sway  and  heave  motion 
of  unit  amplitude  are  given  in  Table  5.  The 
values  given  under  the  heading  "Infinity"  are 
obtained  in  terms  of  the  damping,  B.(,  by 

A<,)  •  £  (B../P)*,  i  -  2,  5 

which  can  ba  obtained  from  the  principle  of 
energy  conservation  and  are  exact  at  y  ■  ”. 

It  appears  that  the  radiating  wave  generated 
by  sway  motion  yields  better  accuracy  than 
that  generated  by  heeve  sot  Ion. 
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IV.  Results  and  Discussion 


The  sample  computation  H  performed  for  the 
condition  of  ze ro  tpood,  l.e.,  0  -  0.  Thus, 
there  Is  no  contribution  from  tha  staady 
translation  to  tha  frea-surface  alavatlon 
which  Is  given  by  Equation  (19). 

To  obtain  tha  ralattva  notion,  tha  absoluta 
vertical  motion  of  tha  body  should  ba  known. 

In  Figures  k  through  t,  tha  notion  amplitudes 
per  unit  Incident -wave  amplitude  ere  shown 
versus  the  nondlnanstonal  frequencies  tu’b/g 
•  2Tb /X  for  u  “  90,  13$  and  l8o  degrees.  Due 
to  the  symmetry  of  the  bod«,  there  Is  no 
motion  at  g  •  w  (head  waves)  for  sway  and 
roll.  To  facilitate  the  understanding  of  the 
measure  of  the  vertical  displacement  at  the 
beam  end  which  is  contributed  by  roll,  the 
roll  amplitude  is  multiplied  by  the  half  bean 
b  in  Figure  6. 


figure  4  —  buy  AmpMude  at  Zara  Spaed  far 
Various  Wave  Heedbtge 
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Since  tha  pontoon  Is  symmetric  fore  and 
aft,  the  aquations  of  motion  given  by  Equation 
(33)  (neglecting  the  surge  motion  end  Its 
coupling  effects)  can  be  reduced  to  the 
following  form; 

(«  ♦  A,,)E,  ♦  ♦  C„5,- 

(H  ♦  AjjIEj  ♦  *t,it  ♦  (A„  -  Hz.K,  ♦ 

♦  -  f1b)  <*» 

(**»  ”  +  *  *»»)£*  *  *»*£» 

♦  ■  F^e) 

Tha  coefficients  appearing  In  Equation  (kl) 
are  obtained  from  Equations  ()t)  and  (3$)  In 
terms  of  tha  two-dimensional  velocity  poten¬ 
tial  0'(y,z;x).  The  (-coordinate  of  the 
center  of  gravity,  («,  was  approximated  as 
that  of  the  canter  of  buoyancy.  The  radius  of 
gyration  for  roll  was  assumed  to  be  0.(b  such 
that  lyy  •  H(0.6b) 1 . 

it  Is  well  known  that  a  prediction  of  roll 
motion  based  on  linear  damping  obtained  by  a 
velocity  potential  at  zero  speed  would  yield 
an  over-predicted  peak  roll  at  Its  resonant 
frequency.  An  Iterative  scheme  to  obtain  a 
convergence  of  roll  motion  using  the  equivalent 
linear  damping  of  the  viscous  damping  could 
have  been  tried  in  order  to  suppress  the  peak 
roll;  however,  to  avoid  the  ambiguity  resulting 
from  the  approximation  of  tha  viscous  damping. 
It  was  decided  not  to  employ  this  scheme. 

The  strip  method  employed  for  the  computa¬ 
tion  of  motion  does  net  take  Into  account  the 
end  effects  of  the  pontoon.  Thus,  the  length¬ 
wise  Integrals  In  evaluating  the  K,  end  I., 
become  simple  Integrals  such  as  *  “ 
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f\t  *  ’ 


\tL 


f\l*2  dx 


m 


where  !*.  tV  sectional  quantity  of  or 
Bu .  Similarly,  tor  the  wave  excited  forces 
Fk(e) ,  the  contribution  from  the  lengthwise 
integral  becomes 


/, 


.JlUcosy  „„ 


2  sln^y  cosy) 
K  cosy 


<*>2a) 


2  xe j  Kxcosy 


dx  -  J 


sln^ycosuj 


L  cos 


■If  -4 


K  cosy 


K*  cos’y 


for  |i  f  r 


(Mb) 


«  0  for  V  «  -y 

For  the  arbitrary  pontoon  geometry  chosen 
in  the  present  sample  calculation,  the  pitch 
and  yaw  motions  for  y  A  a/2  do  not  bear  much 
significance  because  they  are  functions  of  the 
length  of  the  body  which  is  arbitrarily  chosen 
in  the  present  case.  Thus,  the  results  to  be 
illustrated  will  be  limited  to  beam  waves, 
i.e.,  u  ■  a/2,  only.  However,  the  motion 
amplitudes  In  Figures  b  to  4  for  y  A  a/2  are 
presented  to  show  that  for  a  body  «dilch  is 
symmetric  fore  and  aft  as  well  as  port  and 
starboard,  beam  waves  could  produce  larger 
relative  motion  than  the  other  wave  headings. 
It  is  of  interest  to  note  in  Figure  5  that  the 
heave  amplitude  Is  less  than  the  Incident  wave 
amplitude  at  all  wave  lengths  for  both  y  *  13$ 
and  180  degrees.  This  fact  seems  to  reflect 
the  behavior  of  the  function  given  by  Equation 
(b2b)  which  is  the  lengthwise  contribution  of 
the  heave  exciting  force.  One  can  readily  see 
that  the  function  sin  x/x  where  x  ■  cosy  Kl/2 
has  the  maximum  value  of  unity  at  x  ■  0  end 
mono ton I cal ly  decreases  as  x  Increases.  Thus, 
as  either  the  absolute  value  of  cosy  or  K 
increases,  the  wave  excited  heave  force  de¬ 
creases.  The  swey  motion  contributes  to  the 
relative  motion  through  Its  wevemoking  which 
changes  the  free-surfece  elevation  at  the  ship 
sides. 


The  local  waves  generated  by  the  sway, 
heave,  end  roll  motion  ere  presented  In  Figure 
7  as  Aj,  Aj,  and  \  ,  respectively,  where  the 
bar  Indicates  that  the  quantities  ere  non- 
dimensional  I zed  by  Incident  wave  amplitude.  The 
beaemlse  location  of  these  free-surfece  ampli¬ 
tudes  Is  1 .05b  from  the  origin.  As  can  be 
observed  In  the  figure  the  maximum  local  wave 
amplitudes  occur  at  the  same  frequencies  at 
which  the  respective  motion  amplitude  becomes 
maximum.  The  trend  of  the  curves  shown  In 
Figure  7  Is  found  to  be  almost  Independent  of 
the  beaaMlse  locations,  although  the  magni¬ 
tudes  can  change  significantly.  For  the  body 
chosen  In  the  present  case,  the  sway  motion 
appears  to  be  a  better  wevemaker  than  the 
heave  motion  per  unit  amplitude  of  the  Inci¬ 
dent  wave  from  the  beam  directions  for  higher 
frequencies,  I.e.,  Kb  >  0.8. 


Figure  7  -  Maximum  Free  Surface  EleveHen  at 
y  •  ±  1.88b  Duo  to  Motion  In  Durey 
(XjL  Heave  |X,|  and  Hal  I  A,) 


Combination  of  these  motion-generated  waves 
will  be  a  harmonically  oscillating  free  sur¬ 
face.  The  amp  I Itude  of  this  free-surfece  os¬ 
cillation  at  y  •  t I .05b  Is  shown  In  Figure  8 
as  Am  which  is  a  nondlmensional  quantity 
normalised  by  the  incident  wove  amplitude. 

The  free  surface  deformetlon  due  to  the  Inci¬ 
dent  and  diffracted  waves  should  be  added  to 
the  motion-generated  free-surfece  disturbance 
to  obtain  the  actual  free-surfaca  elevation  at 
the  sides  of  a  ship.  The  diffracted  waves  are 
the  reflected  end  the  transmitted  waves.  At 
the  Incident  side  of  the  body,  the  free  sur- 
face  deformation  Is  caused  by  the  combination 
of  the  Incident  end  reflected  waves.  The 
amplitude  of  the  comb  I nod  waves  of  Incidence 
and  reflection  at  the  incident, side  at 
y  ■  1.05b  Is  Indicated  by  In  Figure  8. 

The  transmitted  wave  amplitude  at  the, lee 
side  at  y  •  -1.05b  Is  Indicated  by  V1'.  The 
bar  sign  indicates  a  nondlmensional I  tat  Ion  by 
the  Incident  wave  amplitude.  As  one  can  ob¬ 
serve  from  Figure  8,  the  free  surfece  oscil¬ 
lates  near  the  body  et  the  Incident  side  with 
twice  the  amplitude  of  the  Incident  wave  for 
lb  >  0.4,  I.e.,  the  Incident  wave  length  less 
than  about  10b.  On  the  other  hand,  at  the  lee 
side  the  free  surface  fluctuates  with  less  than 
about  bO  percent  of  the  Incident  wave  ampli¬ 
tudes  for  A  >  0.4. 
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•y  the  principle  of  the  conservation  of 
energy  the  fol loving  relation  should  hold 

k'1  *H  -- 

*  ’ym  m  '  ym  >• 

where  and  ere,  respectively,  the 

amplitude  ratios  of  the  reflected  eed  trans- 
mitted  waves  to  the  Incident  waves.  A  check 
was  made  at  |y|  *  100b  to  confine  the  fore* 
going  relation.  It  was  found  that  the  rela¬ 
tion  holds  true  even  at  distances  as  short  as 
I 00b.  However,  when  the  distribution  of  the 
source  Is  extended  onto  the  waterline  Inside 
of  the  body,  the  square  sum  of  the  foregoing 
equation  becomes  slightly  less  than  unity. 

This  Is  construed  as  an  Indication  that  the 
added  top  deck  may  Introduce  a  slightly 
greater  numerical  error  than  the  case  without 
the  top  deck.  It  Is  uncertain  If  Kg'1! can  be 
greater  than  2  even  at  values  of  y  as  seal  I  as 
1.05b.  Although  no  fine  proof  Is  established, 
It  appears  that  the  values  of  V1'  exceeding 
2  are  due  to  numerical  errors  resulting  from 
the  segment  approximation. 

The  amplitudes  of  the  free  surface  mot  Ion 
generated  by  the  body  motion  and  the  Incident 
and  diffracted  waves  are  shewn  In  figure  9  at 
the  Incident  and  the  lea  sldas  of  the  body. 

The  prominent  hump  and  hollow  trend  of  the 
amplitude  curves  Indicates  the  sensitivity  of 
the  free  surface  movement  with  respect  to  the 
Incident  wave  frequency.  The  drastic  change 
from  the  behavior  of  Individual  wove  amplitudes 
shewn  In  figure  I  Implies  that  the  phase  dlf* 
fe ranees  between  the  motion-generated  waves 
and  the  diffracted  waves  con  vary  from  0 
degree  to  tlM  deg  root. 


Hpwl- 


To  show  the  behavior  of  the  amplitudes  of 
the  waves  as  a  function  of  y,  the  notion- 
generated  wave  and  the  reflected  wave  at 
A  •  0.8  are  chosen  and  are  presented  In  figure 
10.  The  values  shown  ere  those  normal Ixed  by 
CA.  The  solid  curve  Is  for  the  notion-gener¬ 
ated  wave,  and  the  dotted  one  Is  for  the  re¬ 
flected  wove.  A  larger  scale  for  y/b  Is  taken 
near  the  hull  to  examine  the  behavior  more 
closely,  for  kb  •  0.8  the  wave  length  (2)  Is 
7.85b.  One  can  observe  at  about  y  -  12b  *  I .52 
that  the  outgoing-wave  amplitudes  become  Iden¬ 
tical  to  the  fer-fleld  values.  The  local  waves 
appear  to  be  confined  within  one-beam  distance 
from  the  hull.  Generally,  this  trend  Is  con¬ 
firmed  at  other  frequencies.  Hmever,  when  all 
the  waves  are  superposed  with  their  proper 
phases,  the  amplitude  behavior  becomes  much 
Mrs  radical  as  shown  by  the  chained  curve  in 
figure  10.  The  sharpness  of  the  curve  Is  In¬ 
dicative  of  a  possibility  of  breaking  waves 
somewhere  between  y  »  1.5b  and  A. Ob. 


o.S i— j—  i.  a— i— i  l-i _ i _ i  aaej — i — i _ u> 

101112121411  SO  ISO  240  220 
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The  amplitude*  of  the  absolute  (|£.|)  and 
the  relative  (|Tg|)  motion  divided  by  the  Inci¬ 
dent  wave  amplitude  are  shown  In  Figure  II  for 
y  -  1.05b  and  In  Figure  12  for  y  •  -1.05b.  The 
relative  motion  computed  on  the  basis  of  no 
deformation  of  the  Incident  wave  Is  designated 
as  "Old  Rel."  and  the  present  calculation  based 
on  the  deformation  of  the  incident  waves  Is 
designated  as  "New  Rel.'1.  A  large  difference 
between  the  old  and  new  relative  motion  can  be 
observed  in  Figure  II  at  higher  frequencies. 

The  difference  is  almost  a  factor  of  two.  This 
phenomenon  is  due  to  the  fact  that  at  higher 
frequencies  the  absolute  vertical  motion  of 
the  body  diminishes  while  the  reflection 
effect  almost  doubles  the  amplitude  of  the 
free  surface  motion  on  the  side  of  the  body. 

The  old  relative  motion,  however,  assumes  that 
there  is  no  deformation  of  the  Incident  wave; 
hence,  the  relative-motion  anplitude  becomes 
identical  to  the  incident-wave  amplitude.  On 
the  other  hand,  the  new  relative  motion  at  the 
lee  side  diminishes  at  higher  frequencies 
since  both  the  absolute  motion  and  transmitted 
wave  amplitudes  diminish,  while  the  old  rela¬ 
tive  motion  amplitude  becomes  the  same  as  the 
incident  wave  amplitude. 


0  04  0.8  1.2  16  20 

Kb 

Figure  IT  -  Relative  Motion  and  Abeokita  Motion 
at  y  -  1.06b  for  Beam  Wave* 


0  0  4  0.8  1.2  1.6  2.0 


Figure  12  -  Relative  Motion  and  Absolute  Motion 
at  y  »  -1.06b  for  team  Wave* 


V.  Conclusion 

From  the  present  Investigation  It  Is  found 
that  neglecting  the  Incident-wave  deforma¬ 
tion  due  to  the  body  motion  and  diffraction 
could  lead  to  a  significant  error  In  the  pre¬ 
diction  of  the  relative  motion  of  ships.  The 
computation  should  be  extended  to  an  actual 
ship  at  speed  for  various  wave  headings, 
and  the  results  should  be  correlated  with 
available  model  experimental  results  of  rela¬ 
tive  motion. 
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Appendix  I 

Determination  of  the  Velocity  Potentials 
Pp(y,z)  and  pfc(y.z) 

The  boundary  condition  to  be  satisfied  by 
the  harmonic  functions  Pp(y,z)  »"d  PkW.*)  are 
as  follows: 

g  -  op  “0  on  z  •  0  (1-1) 

where  a  ■  wVg  for  P0  and  o  •  w‘/g  for  p  ; 

p„  I  •  f  (1-2) 

'c. 

where  C,  is  the  contour  below  the  calm  water 
surface  of  a  cross  section  of  a  ship  and  f. 
a  complex  function,  is  assumed  to  be  known; 

lim  (g‘  T  jap)  “  0  (I-J) 

y4  ♦»  y 


•  0  is 

The  Green's  function  G(y,z;n,C)  which  Is  in 
the  form 

G  -  In  r  ♦  M(y,z;n,c) 


r  -  l(y-n)1  ♦  (z-c)M* 


($♦$)— 

in  z  <_  0,  I*  given  In  Wehausen  and  Laitone^  as 


G  -  ijin  r  -  tn  r  -  i/ 


♦  J2*ea(m)coso(y-n)|  (1-5) 

where  r'  »  Ky-n)*  ♦  (rtc)1)*.  Using  this 
Green's  function,  w*  can  express  P'  by 

♦  ’(y.t)  •  y*Q(l)GO;y,z)d1  (1-6) 
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* 


where  Q  Is  th«  unknown  to  bo  determined. 

Since  6  represents  e  pulsating  source  at  a 
point  below  the  cal*  Mater  surface,  4  can  be 
regarded  at  an  unknown  source  strength  of 
coaplex  value. 

Following  the  procedure  adopted  by  Frank,16 
we  can  show  that 

V  (y  ,z)  *  £  <7.  f  6(l;y,*)dl  (1-7) 

I -I  7, 

Cl 

where  C(  Is  divided  Into  N  even  ma*er  of 
linear  sage  ante  C ( ,  and  Q|  It  the  constant 
source  strength  on  the  Ith  line  segment  C(. 

Let  us  ass-:-'  >n*t  the  prescribed  function 
f  can  be  if  »  rs-co  an  even  function,  f'*J, 
and  ar;  aii  *••»»  ■>’  .  f  '*'•  of  y.  The  cross 

sectier  bxn*«-«“.-j#lon  bare  Is  assiaaad  to  be 
syersfc  «£•■*»■  •’be  /-axis.  Thus,  wa  can  also 
assuau  <r  in>.  that  ♦'  can  be  divided  Into 
an  eve?-'  ^  ,y’*)  and  an  odd  part  ♦'** .  Than, 

we  can  <,  &■&-.  tfc«st: 

g(tf*  a.'L'j|*y,6(*)(!  ;y,*)dl  (I-l) 

Cl 


£)  q}°/ e(,)dl 

1-1  1  ■'c, 


<!-»> 


where  N  -  M/2  Is  the  masher  of  linear  lagwants 
on  the  right-half  of  C(,  and 


w  )  r 

(t)  >  -  |tn  r  ♦  In  i-!  -  In  r’  ♦  tn  rj 
*  r  •  k(r»c)  .  .  —  . 


»k(-CWiX±3] 

k-a 

♦  Jl*ea(z*C*cosa(y+Ti)] 


/«•  k 

S - dk 

a 


(x-io) 


where  the  upper  end  lower  Signs  correspond 
respectively  to  fil*)  and  and 


r,  -  ((ywi)*  ♦  («-C)*l*. 

r\  ■  ((r«l)*  ♦  <*n)*J*. 

How,  If  we  take  the  no  ratal  derivative  of 
gl*)  and  g'*»  at  the  midpoints  of  C|  and 
designate  theta  by  g|t)and  g|i),  than  due  to 
the  property  of  the  Vurface  llstrlbutlen  of 
source'^,  we  have 


JO -  JO 

*IH  *  rl 


6<*)(|5y,..,)dl 

(I-ll) 


These  two  equations  can  be  expressed  In  alge¬ 
braic  form  es 

,(.),(.)  .  ,(•) 

-  .<*> 

where  A^  and  A***  are  N  by  H  o-et  rices,  the 
elanants  of  which  are  given  by 

Ajf  -/  S*4)  (l;y,.*,)dl  for  (d)  -  (e)  or  (t) 


and  are  the  ooluaet  vectors,  the 
alaawnts  of  which  are  given  by  W**  and  QJ'», 
respectively  (however.  If  k«l,  only  fmlf  of 
the  values  should  be  taken),  and  I'*'  and  »•'' 
are  the  coluam  vectors,  the  elements  of  which 
are  given  by  fj*»  and  f }•> ,  respectively. 

Then,  we  can  obtain  the  unknowns  and 
d1*'  by 


Q<d)  m  JAC<f)|  for  (d)  »  (e)  or  (•) 

and  substitution  of  these  Into  Equations  (I -I) 
and  (1-9)  yields  g (•)  and  g(0. 

For  the  diffraction  potential  (^,  we  have 
f  =  f(«>  ♦  f(,)  -  -  £ 

6e 

hence. 


f(,)  -  Jwt^^t^slnu  sin(Kyslnu) 

*  N,cos(Ky  slnu)] 

f(,)  -  -wt^lH.sInp  cos(Ky  slnp) 

♦  N|sln(Ky  sinp)] 

where  (y,x)  Is  a  point  on  C(  and  Ht  and  H,  are 
respectively  the  y-  and  z -cosmonaut  of  the  unit 
vector  on  C*  pointing  Into  the  body.  For  a 
sywaetrlc  body  Nt  Is  an  odd  function  of  y  and 
Hf  Is  an  even  function. 
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For  the  notion  potential  ,  we  have 
the  following  boundary  condition  on  C(: 


where  M,  *  yM,  -  zNa. 

Appendix  II 
Irregular  Frequencies 

It  is  well  known  that  the  Green's  function 
method”  for  the  velocity  potential  function 
associated  with  an  oscillating  body  In  a  free 
surface  suffers  f ran  the  existence  of  In¬ 
definite  solutions  at  certain  discrete  fre¬ 
quencies.  These  frequencies  are  often  re¬ 
ferred  to  as  "irregular  frequencies."  John” 
showed  that  they  are  equivalent  to  the  eigen¬ 
values  of  the  boundary-value  problem  for  the 
velocity  potential  defined  In  the  Internal 
fluid  domain  bounded  by  the  body  surface  and 
the  interior  line  zj“  0.  That  is,  If  the 
velocity  potential  d  defined  Inside  the  body 
has  nontrivial  solutions  for  the  boundary 
condition 

♦2  *  Ap  -  0  on  CjU  *  0,  -b  <  y  <  b) 

I  n  r  (M-0 

♦  •  0  on  C, 

for  discrete  values  of  A,  then  aii  •  /gX ,  Is 
the  Ith  Irregular  frequency  for  I  •  I,  1,  .  .  . 

If  the  velocity  potential  d  defined  In  the 
fluid  region  R  outside  the  body  Is  expressed 
by 


•J  Q( 1)6(1 ;p)dl, 

cs 


f(p,)  2  dw(p)|  * 


♦j£  «<l>yi;pt) 


"here  f(p,)  Is  prescribed.  The  Fredholm  theo¬ 
rem  states  that  Q  has  a  unique  solution  If  the 
associated  homogeneous  equation,  l.e.,  f  ■  0 
in  Equation  (II-J),  has  a  trlval  solution 
only.  From  John's  uniqueness  proof*  we  know 
that  if  du|r  ■  0,  d  Identically  vanishes  In 
R  and  on  t.'  Since  d  given  by  Equation 
(11*2)  Is  continuous  everywhere  In  i  <  0.  It 
can  be  extended  to .the  Interior  dome In  R'. 
if  we  let  IV)  -j£f  «(»)•( I i»"pR')dl,  ? 


should  have  trivial  solutions  for  the  boundary 
condition  (II-I)  If  d  I*  unique,  lut,  the 
fact  that  for  same  Xj .  J  has  nontrivial  solu¬ 
tions  implies  that  h*nce> 

Q(l)  -  td)^  0  on  C7  l"X  -  Xj.  The  Iden¬ 
tity  Q(l)  ■  dq(l)  can  be  easily  proved.  Frqm 
Green's  theorem  wt  can  show  for  any  point  p 
In  the  Interior  domain  R'  bounded  by  C#  and 
Cf  that 

d(p')  -f  ldN6(p‘)  *  ?6M)dl, 

*'c, 

0-jf  [d||G(p' )  -  dC||]dl ,  (II-4) 

hence,  subtracting  the  two  equations  above,  we 
get 

*<p')  -/*((?„-  d„)6(p*)  -  (d  *  d^Jdl 

mfc  V(p,)dl  (II-5) 


since  d"dadNB0on  C  ,  and  from  Equations 
(II-A)  and  (II-5)  we  find  that  dq(l)  “  Q(D. 
Since  Q(l)  cannot  be  identically  zero  on  C(, 
the  associated  homogeneous  equation  of  (ll*3) 
can  have  nontrivial  solutions  at  X  •  X(. 

Frank*0  has  shown  that  the  added  mass  and 
damping  coefficients  of  cylindrical  bodies  have 
discontinuities  at  certain  discrete  frequencies. 
Several  Investigators  have  shown  either  in  pub- 
II shed"  •*’»*"  or  unpublished  forms  various 
methods  for  alleviating  the  irregular  fre¬ 
quencies.  The  removal  of  the  irregular  fre¬ 
quencies  is  achieved  in  this  work  by  Imposing 
dz  ■  °  or  d'  "  0  on  the  interior  weterline  Cf. 
This  Imposed  boundary  condition  on  Cf  then 
makes  d  =  0  In  R*  since 

Jf  |Vf  |*  ds  -  f  d  ♦„  dl  -  0  (II-6) 

»'  Vcf 

for  J|r-  -  0  end  fplcf  •  0  and,  therefore. 

|V$|  mo  In  R  implies  that  f  *  0. 

With  this  approach  we  can  begin  with  a  new 
definition  of  d*  by 


Q(I)G(I  ;p)dl 


and  eolve  for  Q  from  the  Fredholm  equation 

tip,)  •  •  *<(»*)  ♦  J  <0)«|,(l»t)dl 
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